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This contribution - a continuation of the work in refs.[Q]-Q - reports on recent developments in the calculation 
and understanding of 4-quark energies generated using lattice Monte Carlo techniques. 



1. THE CALCULATION OF 4-QUARK 
ENERGIES IN A TETRAHEDRAL 
GEOMETRY USING MONTE CARLO 
SIMULATION 

In refs.^,^ - the binding energies of four 
quarks in various geometries (e.g. at the cor- 
ners of a rectangle or on a straight line) were 
calculated in static-quenched-SU(2) on a lattice 
16^ X 32 with 13 = 2.4. This showed that the 4- 
quark binding energies were greatest when two of 
the three possible partitions of the four quarks - 
into two colour-singlet pairs - were degenerate in 
energy e.g. as in a square geometry. The main 
reasons for studying other geometries were: 1) to 
investigate further this point and 2) to give a rep- 
resentative set of geometries illustrating the dif- 
ferent possibilities for arranging four quarks. The 
latter would then serve as a very stringent test 
of models attempting to describe four interacting 
quarks - provided such models involve approxi- 
mations analogous to those contained in the cor- 
responding lattice calculation e.g. static quarks 
and SU(2). 

Interest in the tetrahedral geometry arises, 
since in this case all three partitions are degen- 
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erate in energy. For completeness, a series of re- 
lated geometries, where the four quarks have the 
coordinates (0, 0, 0), (r, 0, d), (0, d, d) and (r, d, 0), 
are now studied. Preliminary results are given in 
table 1. 



Table 1 

The ground and first excited state energies (i?o,i)- 
The results for r = are from ref. [4] for the 1x1- 
square using 2 partitions and those for r ^ d for 
the Tetrahedron. All energies are given in terms 
of lattice units dX (3 — 2.4, which corresponds to 
an inverse lattice spacing « 1.7 GeV. 
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-0.0696(2) 
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-0.016(2) 
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-0.043(3) 
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0.085(2) 
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0.147(2) 
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(3,3) 


(3,4) 
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-0.041(1) 


-0.028(1) 


-0.010(2) 
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0.048(2) 


-0.028(1) 


0.105(5) 
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(4,3) 


(4,4) 


(4,5) 


So 


-0.039(2) 


-0.036(1) 


-0.012(2) 


El 


0.03(1) 


-0.036(1) 


0.089(4) 



Here two points are noteworthy: 

1) For the tetrahedral case (i.e. r = d), the 
ground state energy {Eq) and the first excited 
state energy (Ei) are degenerate. As will be seen 
in section 3) below, this degeneracy could have a 
natural explanation using a simple model. 

2) As r becomes small, the results match on to 
the earlier ones in ref. for r = - the d x d 
square geometry, where only the two degenerate 
quark partitions were employed to give the bind- 
ing energies. For example, with {d,r)= (2,1) and 
(3,2) the values of Eq are the same when using 
either the two degenerate partitions or the com- 
plete set of three. However, the results for Ei 
with two partitions get less accurate as r increases 
giving 0.086(2) and 0.053(4) for the above values 
of {d,r). This gives added confidence that, for 
squares, the earlier use of two partitions was cer- 
tainly adequate. 



2. EXTRACTION OF BINDING ENER- 
GIES USING THE FULL CORRELA- 
TION MATRIX 

In refs.H,^ the energies were extracted by first 
solving for the eigenvalues of 

WLo'l ^ XPw^-'al. (1) 

Here, for euclidean time differences T, the 
are the values of Wilson loops with m, n = 1,2,3 
for the 2-quark energies (i.e. 3 fuzzing levels 
12,16,20) and m, n = 1, 2 or 1, 2, 3 for the 4-quark 
energies depending on the number of partitions 

(T) 

used - all with fuzzing level 20. Knowing A^^ ' , the 
2-and 4-quark energies are given by the plateau 
as Ei = — lim7-_^oo In ^A'"^''^ . In practice, since 
it is the 4-quark binding energy that is of prime 
interest, the plateau is extracted directly for the 
difference of the 4- and 2-quark energies. These 
differences are quoted in table 1. The errors are 
then estimated by bootstrapping this procedure. 
As discussed in ref.|5|, this technique does not 
make full use of the available Monte Carlo data, 
since it disregards correlations between different 
values of r. A more complete fit minimizes the 



following function 

DiT'P')C-^D{TP), (2) 

T,T',P,P' 

where D{TP) = [F{TP) - W{TP)] . Here W is 
the average of the Wilson loop data and F is a 
function to be fitted to this data. Both of these 
are dependent on T and also P - an index denot- 
ing the fuzzings or partitions involved i.e. correla- 
tions among the data in both T and in P are now 
included. C{T'T, P'P) is the full correlation ma- 
trix among the measured Wilson loop data sam- 
ples. Since inverting a large matrix with limited 
data samples is unstable, we need to approximate 
C by some model e.g. 

CiT'T, P'P) = A{P'P) exp(-a|T' - r|) 
would be one particularly simple choice. Usually 
F has a form such as a,i{P) exp{—miT) , where 
the nii are then interpreted as the energies (Ei) of 
the system. Such an analysis is now being carried 
out. Preliminary results in table 2 indicate that 
the earlier energies and their error estimates were 
indeed correct. 

Table 2 

Comparison between the energies Ei{C) ex- 
tracted using the full correlation matrix and the 
earlier results Ei(P) with the plateau method. 
These results are for quarks on the corners of rect- 
angles with sides (d,r). 



(d,r) 


(2,2) 


(2,3) 


(2,4) 


EoiC) 
EoiP) 


-0.0590(3) 
-0.0585(4) 


-0.005(1) 
-0.006(1) 


0.0006(7) 
-0.0004(4) 


EiiC) 
Ei{P) 


0.141(1) 
0.143(1) 


0.325(2) 
0.324(1) 


0.514(1) 
0.512(1) 



3. A MODEL FOR UNDERSTANDING 
4-QUARK ENERGIES 

In refs.||]-^] a model for attempting to under- 
stand the above lattice results was introduced. 
This amounts to comparing the lattice energies 
with the eigenvalues [Ei{f)] from 

[V(/)-i?,(/)N(/)]*, = 0, (3) 



with 
N(/) = 

V(/) = 



1 //2 
//2 1 

Vl3 + V2i fV, 



BA 



AB 

Via + W23 



Here Vab Vba ^ 

- (Ul3 + V2i + Wl4 + W23 



Vl2 - W34) 



(4) 



(5) 



as expected from a t^.Tj colour- vector potential. 
For squares of a given size c? x d, one value of f{d) 
gives a reasonable fit to both Eo{d) and Ei{d). 
However, for the tetrahedral geometry, since all 
the Vij are equal, the two eigenvalues of eq.(|^) 
are degenerate at zero binding energy. To push 
this degeneracy to some other energy could be 
explained by including a colour-scalar interquark 
potential in addition to the colour-vector poten- 
tial used to derive eq.(5). 

Further support for the above model comes 
from ref . , where it is shown explicitly that per- 
turbation theory upto fourth order in the quark- 
gluon coupling does indeed correspond to the 
/ = 1 hmit. 

Of course, the above replacement of the i?o,i 
by / is not particularly useful unless / - a func- 
tion of all the four quark coordinates - can be 
parametrized in some convenient and simple man- 
ner. So far, the most successful appear to be the 
forms 



/ = /cexp[-a6oA - ^^/VaP], 



(6) 



where &o is the string tension dictated by the 
2-quark potential, A is the minimal area of the 
surface bounded by the straight lines connecting 
the quarks, P is a perimeter connecting the four 
quarks and a, 7 are phenomelogical constants ad- 
justed to fit the lattice data. In this way, the orig- 
inal rapid variations of £'0,1 with geometry can 
be first described in terms of a smoothly varying 
function /, which can then be expressed in terms 
of the constants 0,7. In practice, two specific 
parametrizations have been attempted. The first 
assumes /c = 1 - consistent with the expec- 
tation that / ^ 1 as A and P become small i.e. 
in the perturbative limit. In addition, guided by 



arguments based on the lattice cubic symmetry, 
P is dependent on the underlying lattice - not 
just the quark positions - and so, being a lat- 
tice artefact, must be dropped in the continuum 
limit. In the second parametrization, P is defined 
simply in terms of the quark coordinates by the 
straight lines joining the quarks. However, /c is 
taken to be (3 dependent. A preliminary fit gives 
/c(/3 = 2.4) = 0.88(2) and Ufi = 2.5) = 0.94(2) 
i.e. apparently /c — > 1 as /3 increases - again as 
expected in the perturbative limit. The contin- 
uum form of / then is taken to have /c = 1, with 
A and P defined purely by the four quark posi- 
tions. At present, it is not clear which of these 
two parametrizations is to be preferred. 
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